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Abstract 

We calculate the two-step contribution to (p,p') and (p, n) reactions at intermediate energy. We 
describe the motion of the incident nucleon with a plane wave and compare the contribution from 
two-step processes with that from one-step processes. To describe the two-step processes, we extend 
the response functions to nondiagonal forms with respect to the momentum transfer q. 

We performed a numerical calculation for the cross sections of the 12 C, 40 Ca(p,p') scattering and 
^ the spin longitudinal and spin transverse cross sections of the 12 C, 40 Ca(p, n) reactions at 346 MeV 

and 494 MeV. We found that the two-step contribution is appreciable in comparison with the one-step 
■ processes in the higher-energy transfer region for the spin longitudinal and the spin transverse (p, n) 

reactions. We also found that the two-step processes give larger contributions to the spin transverse 
(p, n) reaction than to the spin longitudinal reaction. This finding is very encouraging to interpret 
the discrepancy between the DWIA calculation and the experimental results of the spin longitudinal 
O ' and the spin transverse cross sections. 

> 

X. 

contribution to the cross sections. When the excitation energy is low, it is probable that the projectile 
nucleon collides with the nucleons in the target only once, but as the excitation energy becomes higher, the 
pre-equilibrium processes, in which the projectile collides with the nucleons several times, arc considered 
to become more effective. 

Actually, it has been reported that the single-step process calculation underestimates the scattering 
cross sections, especially in the large scattering angle region, Jy, ^| where the two-step and further multi- 
step scattering are found to have a large effect. 

The multi-step direct reaction (MSDR) was actively studied at the end of the 1970s. Feshbach, Ker- 
man, and Koonin (FKK) j^, |j developed the framework of multi-step reaction theory. Tamura, Udagawa, 
and Lenske (TUL),Q| who pointed out a problem in FKK theory, replaced the sum over the excited nu- 
clear eigenstates with that over 1-particle-l-hole (lp-lh) states, introducing a weight function. Another 
type of formalism was presented by Smith and Wambach, j^] who applied the Glauber approximation to 
the motion of the projectile and analyzed the forward angle scattering. Kawai et al. ^, ^, || presented 
a multi-step formalism with a local-density (or Thomas-Fermi) approximation for the target and the 
semiclassical distorted wave (SCDW) approximation for the projectile and the ejectile. They calculated 
the cross sections of 58 Ni, 90 Zr(p,px) and (p,nx) reactions at incident energies lower than 200 MeV and 
determined the scattering angle dependence of the one- and the two-step processes. They extended their 
calculation to three-step processes [|| and obtained comparable results with one- and two-step processes 



1 Introduction 

In nucleon induced high energy reactions, single-step processes have often been considered as the main 
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at large angles. De Pace et al. (TOH also used the Glauber approximation for the analysis of (p, n) reactions 
and took account of the spin transfer. 

Although multi-step reactions have been studied by various groups, as mentioned above, those with 
spin dependence have been rarely investigated. In the spin polarized cross sections of the (p, n) reactions, 
a large discrepancy between the distorted wave impulse approximation (DWIA) calculations and the 
experimental results has been reported. 

In the (p, n) reactions, the spin longitudinal response function Rt and the spin transverse response 
function Rt are extracted from the polarization transfer coefficients, jy], |l2|, O, |l4| With these observ- 
ables, the ratio Rl/Rt is found to be less than 1. 

It has been predicted with the random phase approximation (RPA), however, | ji"5| that Rl is enhanced 
and the peak of its energy spectrum is shifted downwards, while Rt is quenched and its peak is shifted 
upwards. The ratio Rl/Rt becomes greater than 1 theoretically, and this contradicts the experimental 
results. 

The spin longitudinal cross section ID q is roughly reproduced by DWIA with the RPA correlation 
, in the lower excitation energy region, but in the spin transverse cross section ID p , the estimation 
amounts to only about half of the experimental result in the whole excitation energy region (see Figs. 1,2). 
In the region where the cross sections are underestimated, the multi-step processes are expected to have 
some contribution. 

For this reason, we calculate the two-step spin longitudinal and the spin transverse cross sections of 
the (p, n) reactions. These reactions were already studied by De Pace et al., but they applied the Glauber 
approximation, and they assumed that the first step momentum transfer and the second step momentum 
transfer are nearly parallel to the total momentum transfer. They also assumed that the spin flip occurs 
only in the first or second step and that in the other step the spin-isospin scalar transition takes place. 
Further, they neglected the interference between processes with different shares of the momentum transfer 
in the first and the second steps. 

In this paper, we develop a formalism for two-step reactions within the framework of the plane-wave 
approximation in order to see the relative contributions from one- and two-step processes. We believe 
that the path lengths of the incident particles in the target nucleus in the one-step and the two-step 
processes are nearly equal, and therefore the effects of absorption in the one-step processes and in the 
two-step processes become nearly equal. We treat the spin degree of freedom carefully, taking account 
of the difference between the direction of the total momentum transfer and that in each step. With this 
formalism, we calculate the unpolarized cross section of the (p, p') scattering and the spin longitudinal 
and spin transverse cross sections of the (p, n) reactions at intermediate energy. We investigate the spin 
dependence of the two-step contribution and partly explain the contradiction in the Rl/Rt problem. 
In §|2| we present our formalism for one- and two-step processes in detail and study some relations with 
TUL's formalism. In we describe the cross sections with spin degrees of freedom. In §0 we explain our 
method for the numerical calculation. In §|| the results and discussion are given, and in §H we summarize 
this paper and give the conclusion. 



2 Formalism 

We consider nucleon induced inelastic scattering and charge exchange reactions in which the target is 
excited to the continuum region. The total Hamiltonian H is 

H = H a + V, (2.1) 

where Hq is the unperturbed Hamiltonian and V is the residual interaction between the projectile and 
the target. They are specified in the projectile-target center of mass (cm.) system as 

H a = H p + U + H Tl (2.2) 
H p = v/MHF, H t = M T + H%* + ^- , (2.3) 

A 

V = J2m~U =V-U, (2.4) 

1=1 
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= H^ en + V T , (2.5) 



.1 



^2 



i— 1 L J i<j i— 1 

where H p and -f/r are the projectile and the target Hamiltonian, respectively, and v is the effective 
interaction between particles. The subscript represents the projectile and i designates the i-ih nucleon 
in the target. U is the mean field that the target nucleus creates for the projectile. M and Mt are the 
nucleon mass and the ground state target mass, respectively, p is the projectile momentum operator in 
the projectile-target cm. frame, and pi is the momentum operator of the i-th nucleon in the intrinsic 
frame of the target. The XJ[ are the mean fields for the nucleons in the target. 
The T-matrix is given in the 'solved' Lippman-Schwinger form 

T = V + V E+ \ H V, (2.7) 



where E + = E + ie. Using (2.1) and expanding l/(E + — H) in V, one obtains 



with 



Then, the T-matrix elements become 



T = J2V(G Vy-\ (2.8) 

3=1 



T^ = Y, T no\ (2-10) 

3=1 



with the j-th order T-matrices 

= ( X «|($ n |y(G'oFr 1 |$o)|x i (+3 ), (2.11) 

where x[ + ^ an( l Xi "* are the distorted waves of the initial and final channel, and the eigenstates of the 
target nucleus |$„) obey the equation 

H^ n )=E^ n ). (2.12) 



Here n = denotes the ground state, and we set E^ 11 ' = 0. The j-th term of (2.10) describes the j-step 
processes. 

Here we neglect the interference terms with j ^ j' . The different order T-matrices mainly excite 
different number particle-hole states. Then, the final state of is different from that of '. Hence 
the contribution from the interference terms can be neglected. Then the double-differential cross section 
becomes 

^-KY X« (2 13) 

3=1 3=1 

^"(Sl' X (j) =52\T^\ 2 6(u-(Ei-m (2-14) 

where lu (= — ef) is the energy transfer, fc; and kf are the projectile and ejectile asymptotic wave 
number, /Zj (/if) is the reduced energy of the projectile (ejectile), 

Mi = ^' Mf = ^T^' (2 ' 15) 

and ei (ef) and E[ (E{) are the projectile (ejectile) and the target (residual) nucleus energies, 



a = y M 2 + kf, e f = y M 2 + kf, (2.16) 
Ei = M T + ^, Ek-Mv+EP + S;. (2-17) 
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2.1 One-step processes 



We are mainly interested in the relative importance of the one- and the two-step processes. Therefore we 
adopt a i 
becomes 



adopt a plane- wave approximation for the distorted waves Xi an d Xf Then the one-step T- matrix 



T«(k f ,ki) = (kfK* n |^|*o)|ki) = (kf|($ n |V|$o)|ki) (2.18) 
for n^O. The effective interaction V between the projectile and the target can be written as 

A 

V(r ) = w oi(r - r t ) = / v(v - v)p(v)dv, (2.19) 

i=l J 

where p(v) = X^<=i <-K r — r 0- For the moment we suppress the spin-isospin dependence of the interaction 
v to make the argument simple. From §3.1 we explicitly introduce the spin dependence. 

In the impulse approximation the effective interaction v(tq — r) between the projectile and the nucleon 
in the target is approximated by 

»(ro - r) = J *NN(pK p - (ro ~ r) ^, (2.20) 

where £nn (p) is the nucleon-nucleon transition matrix (NN i-matrix) for the momentum transfer p. Then 

(kf|V(r )|ki) - t NN (k f - ki)p(k f - ki). (2.21) 
Defining the momentum transfer as q = kf — kj, one finds 

T$(kf,k0 = tN N (q)<*n|p(q)|*o>, (2-22) 

where p(q) is the density operator in the momentum space. This is the T-matrix in the tp approximation, 
which is a simple version in various impulse approximations. The NN t-matrix is, in fact, expressed as 
a function of q, Q(= k + kf ), and the incident energy in the laboratory frame -Kiabj but we suppressed 
the Q dependence, and we do not write Ki a b explicitly. 
Substituting (2.22) into ( [2.14 ), one obtains as 

X W = |*NN(q)| 2 £ |($„|p(q)|$o)| 2 ^ - {Et - Ek)), (2.23) 

n=iO 

= ^-|tNN(q)| 2 Ji(q,a; int ), (2.24) 
where i?(q, w mt ) is the response function of the density fluctuation 

R(q, c int ) = £<$ |p t (q)|$n)($n|p(q)|$o)5(^ ,lt - E^). (2.25) 

The quantity M R (= M T + E 1 ^) is the mass of the residual nucleus, and s = (e ; + Ei) 2 . The response 
function defined in J2.25 ) is the strength per unit energy with respect to the energy transfer in the 
intrinsic frame ui lnt , while the double-differential cross section is the cross section per unit energy with 
respect to the energy transfer in the projectile-target cm. system ui. The factor \du> mt / du>\ = t/s/Mr is 
the variable transformation coefficient derived from the equation Jl8| 

„ - i±4^M. ,2.26, 
2vs 

Thus the one-step double-differential cross section can be written as a product of the square of the NN 
t-matrix and the response function. 

Since the exact n uclear states |$ n ) are complicated, we cannot calculate the response functions by 
using the expression ( 2.25 ). To obtain a calculable expression for the response functions, we adopt two 



4 



approximations of TUL.[Q They assumed that |$q) is the O-particle-0-hole (Op-Oh) state and divided the 
excited state into lp-lh states |$b) and other more complicated states and wrote it as 



With the equation 



l^n) = $>£ ) |*b) 71 ? 0. 



($ n |p(q)|$o)=E a s ) *(^|p(q)|$o), 



(2.27) 



(2.28) 



which implies that the density operator p(q) can excite the target only to lp-lh states, ( 2.25 ) becomes 

(n) (n)* , 



n=i() B'B 



Then we use the TUL statistical approximation 



E (n)* (n) r/ int pint\ _ r / ,int\ 

njtO 



(2.29) 



(2.30) 



by as suming that the are random with respect to n for a fixed B. For B = B' , the left-hand side of 
( 2.30 ) becomes CB(a/ nt ), the probability per unit energy that the state \&b) is mixed in at the excitation 
energy u> mt . Then we obtain the relation 



i?(q, w int ) ~ ^ Ci3 (a; int )($o|p t (q)|*B)($B|p(q)|^o 



(2.31) 



B 



This is calculable with a reasonable assumption for ce(w lnt )- 
2.2 Two-step processes 

The two-step T-matrix in the plane wave approximation can be written as 

T^(k f ,k ; ) =^{k t \{$ n \V\*n>)G n >{*n>\V\*o)\*i), 



where the propagator G„> is 



1 



E+~H 



We rewrite ( 2.32 ) as 

Til\k u k i ) = J2 J (2^)3 (2^ (kf|($n|V|$ n >)|l4)(km|gn>|k m )(k m K$ n ,|V|$o)|ki; 



(2.32) 



(2.33) 



(2.34) 



Since we are interested in high incident energy reactions, we can assume that the mean field U has a 
small effect on the particle. We assume that the momentum of the propagating particle does not vary 
during the sequential collisions, and thus 



U = U (k m )S(k m - k'J. 
The Green's function becomes diagonal with respect to k m , so we can write the T-matrix as 

T$(k f ,ki)= ]T / ^(k f |($„|V|^)|k m )G„'(k m )<k m |<<i>„,|V|<i> )|k i ), 



where 



G„'(k m ) = 



E+ U (k m ) - ^/M^ + kl - ( M T + Eff + ^J- 



(2.35) 



(2.36) 



(2.37) 
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Now we use the two TUL assumptions. The first one is 

(<MV|*o) = Y, a c h (*cWo). (2.38) 

c 

The other is the 'spectator assumption', 

{<£ n \V\<S> n ,)=J24i*(<S> B \V\<P ), (2.39) 

B 

which implies that a collision always occurs by creating a new particle-hole pair and that the particle-hole 
pairs whic h are present in |$ n ') play only the role of spectators. Substituting ( |2.38 ), (2.39) and (2.21) 



into ( 2.36 ) we obtain 

r^(k f ,k ; )= J ^3«^NN(q-qi)(^B|p(q-qi)|$o) 



BC 



xG„,(ki + qi )a£ ; •tN N (qi)(*o|p(qi)|*o), (2-40) 



where qi = k m — kj . 

Next we proceed to calculate of ( |2.14| ). It becomes 



X (2) = Vly §{uj int _ S int } y y y y f pi f p 



^q'l (")* («) («')* («") 



3"Bn' "B'n'"*C "C" 

x < NN (q - qi)4 N (q - qi)<*olp f (q - qi)l*B'>(*s|p(q - qi)|$o) 

x G* „ (k + qi)G n , (k + qi )i NN ( qi )4 N (qi ) <*„ |p f (qi ) | $C ) <*c |p(qi) I *o> • 



Substituting the TUL assumption 



which is analogous to ( [2.30| ), we obtain 



x (2) = V± y y [ d <li [ d q'i 



(2.41) 



J] ^(^ int - 4 nt )4"i*4t" - WWcs* (2.42) 



x CBn'{^ mt )a { c ) *ac' )i NN(q- qi)4 N (q- qi) 

x ($ |p f (q - qi)|*B)<*s|/»(q - qi)|*o)G«'(ki + qi)G;,(k; + qi) 

x tNN(qi)4 N (qi)($o|p t (qi)| < i>C')( < i > c|p(qi)| < i>o). (2.43) 

With the identity J djjj\8{u\ — (E m — E\)) — 1, X^ 2 ) becomes 

x m = Vly y j^Jp^jpi 

x 5{u x - (E m - E i ))c w (o; int )a^' ) *a^ ) 

x tNN(q-qi)4N(q-qi)(*o|p t (q-qi)|3 > B)($s|p(q-qi)l$o) 

x G(ki + qi , Wi)G* (k; + qi , Wl )i NN ( qi )4 N (qi ) ($„ |pt (qj ($ c |p(qj )|* ), 

(2.44) 

where i? m is the target energy in the intermediate state, 

E m = M T + E^ + ^-. (2.45) 



6 



Here we have eliminated the n' dependence of G n ' (k; + qi ) as 

1 

E+ - C/ (k m ) - ^AP + - (m t + A- + w x ) 



G(ki + qj,wi) = t ^ -. (2.46) 



Rewriting the delta function with lo\ , ( 2.44 ) becomes 

y(2) _ S /" j, , /" foil f jgh r/, ,int pint-i„ A ,int\_(n')* (n') 

K „/^n unci •> J \ i J \ i 



R n'jtOBCC 

^NNV 



x i NN (q - qi )4n» (q- qi)(*o|p t (q- qi)|$s)(*s|p(q- qi)|*o) 



(2.47) 



x G(k i + q 1 , Wl )G*(k i + q 1 , Wl )t NN (q 1 )4 N (q 1 )(ci>o|p t (q 1 )|$ C 0( ( i > c|p(qi)|$o). 
Using the approximation 

^(^^-csC^-wf), (2.48) 

which results from the assumption ( 2.39] ) that the particle-hole pairs which are present in the intermediate 
state have nothing to do with creating a new particle-hole pair, we apply the TUL approximation (2.30) 
as 

x <NN(q- qi)4j N (q- qi) 

x J2 - O^ol/^q - qi)|$ B ) <$ B |p(q - qi)|$o) 

B 

x G(k i + q 1 ;wi)G*(k i + q' 1 ;wi) 

x <NN(qi)4 N (q' 1 )^ccK nt )(<i>o|p t (q'i)|<i>c)<<i>c|p(qi)|$o). (2.49) 

c 

Consequently, we obtain 

dqi f dqi 



Ml J 'J (2tt) 3 J (2tt) 3 
x t NN (q - qi )4 N (q - q'i )i?(q - qi , q - q[ ; w int - u4 nt ) 

x G(k i + q 1 ; Wl )G*(k i + q' 1 ; Wl )t NN (q 1 )4 N (q' 1 )i?(q 1 ,q' 1 ;< t ), (2.50) 
where i?(q, q'; w mt ) is expressed as 

i?(q,qV nt ) =^CA( W int )($o|p t (q')l < fA)(^|p(q)|<5o>. (2.51) 

Here we have extended the respons e fun ction to a nondiagonal form with respect to the momentum 
transfer q. The last three factors in ( 2.50| ), which is the product of two NN i-matrices and the response 
function, represent the first collision, and with the next product of the Green's functions, the particle 
in the intermediate state propagates in the target nucleus. Then the second collision occurs with the 
remaining three factors. 



3 Formalism with spin 

We are concerned with the unpolarized cross section of the (p,p r ) scatterings and the spin longitudinal 
and spin transverse cross sections of the (p, n) reactions. Hence we must treat the spin of the nucleon 
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explicitly. First we give the form of the unpolarized cross section and define [ fL7| the spin longitudinal 
and spin transverse cross sections. We introduce the unit vectors 

q „ ki x k f „ „ 

q = R' n = ^7kH' p = qxa (3 ' 1} 

The scattering T-matrix is generally written as 

T(k f , ki) = T + T„o- • n + f q a a ■ q + f p a ■ p, (3.2) 

where do is the spin operator of the incident nucleon and the Tj are the operators of the target. In this 
section we suppress the isospin degree of freedom. 

The unpolarized differential cross section / is given by 

/ = — — - -KTrTr' \TT^}, (3.3) 

2(2J T + 1) L J V ; 

with the target angular momentum Jt- The symbol Tr represents the trace of the spin substates of the 
scattering nucleon and Tr' denotes the summation of all allowed initial and final states of the target, 

Tr'[TTt] = '£(§ \Ti\$ n )(§ n \T\§ )8(u - {E t - Ek)). (3.4) 

0,n 

The observables Di introduced by Bleszynski et al. |l9| arc 

ID 
IDn 
ID q 

IDr, 



2J T 4 


- 1 


1 




2J T -f 


- 1 


1 




2 J T 4 


- 1 


1 




2J T 4 


- 1 



KTr'[f f$], 


(3.5) 


KTr'[f n fi}, 


(3.6) 


KTr'[f q fl), 


(3.7) 


KTr'[f p f£}. 


(3.8) 



We refer to ID q and ID p as the 'spin longitudinal' and 'spin transverse' cross sections respectively. 
PI [l3| In the following, we only consider the case Jt = 0, because we assume that 4>o is the Op-Oh state. 

3.1 One-step processes 

The NN i-matrices can be decomposed generally as 

*NN(q) = }^<ro fi (rip,t^(q), (3.9) 

up, 

with fJ.,ft — u,q,n,p, where o j u = Ij,aj^ = <Jj ■ /t (fi ^ m), (j = 0, i). The operator at denotes the spin 
operator of the nucleon in the target. With this equation, ( 2.24 ) can be rewritten as 

4> = ^E^v^^^^^WCq): (3-ii) 



where 



RM^^) =^c b (c^){$o\pI,(<C\*b){*b\ppW\*o), (3-12) 

B 
A 

p-^^^e-^aip. (3.13) 



i=l 
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The unpolarized, the spin longitudinal and the spin transverse cross sections, /, ID q and ID pi can 
be rewritten as 

I = KJ2K% ID q = KX$, ID P = KX%, (3.14) 

respectively. 

3.2 Two-step processes 

In the spin dependent two-step formalism, in (2.5C) becomes 

J V (2tt) 3 J (2tt) 3 
x S £ X S 2 m>i ct °K ct o^^Om2^o m1) (3-15) 

R M2M2 Ml A 1 ! 

x G(k i + qi;wi)G*(k i + q' 1 ;wi) 

x ^(q-qi.q-qijw^-o^JiZ^Cq^qijwS"*), (3.16) 
where Hx,fii = u,qi,n\,px, and 1^2,1^2 = u, q%, n2,P2, and so on. Here we defined the unit vectors 

q-qi , kfXq2 

q2 = I r, n 2 = —, p 2 = q 2 x n 2 . (3.17) 

|q-qi| |k f xq 2 | 

We also extended the response function to nondiagonal form with respect to the spin direction as 

JV^q';^) = ^c^K lt )($ |pt, (q ' ) |$ A )(< i > A | Pfi ( q )|$ ). (3.I8) 

A 

The unpolarized, the spin longitudinal, and the spin transverse cross sections become 

I = K l -^X^=K [dm iPk, 
2 J J (2tt) 3 J (2tt) 3 

x Yl 12 X ^ 2 i 2 ^ 1 2 Tr ( aa ^< T W2 a ^ <To ^> ( 3 - 19 ) 

dqi /" dqi 



ID =K duJi 



(2tt) 3 7 (2tt) 3 

H l] X S2M>i5^( CT °' 1 i CTo ^^)5^( (T 8 CT o M2 o-o M x). (3-20) 



X 

M2M 2 Ml Mi 



ID P = K [duifp^f-^ 
p J l J (2^) 3 i (2^) 3 

51 51 x S 2 KMi2 Tr(fTo ^' Jo ^ f7p) 2 Tr(CTpCTOp2CTOMi) ' (3 ' 21) 



X 

M2M 2 Ml Mi 



respectively. 



4 Practical method of calculation 

4.1 Various processes of two-step reactions 

We calculate (p,p') inelastic scattering and (p, n) charge exchange reactions at 346 and 494 MeV. In the 
two-step (p, n) reaction, the (p, n) charge exchange collision occurs in the first step or in the second step 



9 



(see Fig. 3). If the (p,n) charge exchange collision occurs in the first step, the second step has to be a 
(n, n') collision, and if the first step is a (p,p') collision, then the second step has to be a (p, n) collision. 

Moreover, there are two cases in non-charge-exchange {p,p') and (n, n') collisions, because the nucleon 
in the target nucleus that participates in the collisions can be a proton or a neutron. Consequently, we 
must consider four cases in the two-step {p 1 n) reaction. 

For two-step (p,p') scattering, we must consider five cases (see Fig. 4). If the (p, n) charge exchange 
collision occurs in the first step, then the second step has to be a (n,p) collision. In this case, the nucleon 
in the target nucleus that takes part in the first step collision has to be a neutron and that in the second 
step collision has to be a proton. This is the only case in which charge exchange collisions occur in the 
two-step (p,p') scattering. 

However, if no charge exchange collision occurs, the nucleons struck in the first step and the second 
step can be a proton or a neutron. Hence we can consider four cases in the two-step (p,p f ) scattering 
that include no charge exchange collisions. 



4.2 Green's function 

Our aim is to compare the two-step contribution with the one-step contribution. We adopted the plane 
wave approximation for the motion of the projectile and the ejectile. Thus in the one-step processes 
no absorption was included. We believe that the path lengths of the incident particles in the target 
nucleus in the one-step and in the two-step processes are nearly equal, and the effects of absorption 
become nealy equal in these processes. Therefore we should also remove the effect of absorption from 
the Green's function. Then the effect of absorption in the two-step processes become zero, as in the 
one-step processes. We use the relativistic expression for the Green's function and apply the on-energy 
shell approximation as 



G(k m ) ~ -iirS \E-U - y/M 2 + k? n - y/(M T + <') 2 + kg, J . (4.1) 

We take a representative value for k m in Uq and set U$ — —3.5 MeV[^0| and 5.5 MeVpH for = 346 
MeV and 494 MeV, respectively. 



Substituting (4.1) into (3.16), we obtain, for instance, from (3.19) 



I = K duj! 



J^d qi / d(cOB0i) / #1 / 7^3^ / d(cOS0i) / d^ 

( Z7T ) J-i Jo Jo ( Z7r ) J-i Jo 

x iwi W o Tr ^ CTo ^ C7o ' i2(7o ' u ^' (42) 



X 



where 



x (2) y 



x i K ^, q'O*^ ( K U> q - qi) W 2 ( K L^ q - qi)Wi ( K ^, qi 

X S \E- Uo - v/M 2 + (ki + qi) 2 - y/(M T + cj[ nt ) 2 + (ki + qi) 2 ^) 
x 6 (e- U - .^Af 2 + (k; + q'J 2 - yj (M T + ^l nt ) 2 + (k; + q'J 2 



x ^(q-q^q-q'ii^-OJJ^jfqi.q'^r'). (4.3) 

We set z-axis parallel to the direction of k ; . With the delta functions, we can carry out the integration 
of cos#i and cos6^ analytically. 

Other integrations are carried out numerically. The mesh sizes were chosen as Aw™* = 10.0 MeV, 
Agi nt = 0.2 (1/fm) and A^i = 15.0°, where q[ nt = {(A - l)/A}q x , with which we can obtain sufficiently 
accurate results. The integration range of qf 11 is restricted up to 4.4 (1/fm), beyond which the response 
functions are negligibly small. 
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Here we have written the energy dependence of the NN ^-matrices explicitly and assumed that the 
projectile has kinetic energy K\ a i, in the first collision and K[ &h in the second collision, where 

KLb = Ki ah - wii ab , wuab = 2^ , t = u{-q{. (4.4) 

4.3 NN t-matrices 

We use the method of Love and Franey |2^| to calculate the NN t-matrix. The derivations of the formulae 
are given in detail in the Appendix. We use the parameter values of Ref. [p3| and treat the exchange 
terms by replacing Q cm with h u as is suggested in Ref. 



A cm 



The scattering amplitude /NN^iab, 9cm) is written as [ 

/NN<Xlab, 9cm) = (Aq + A\(tq ■ Tj)) + (B + Bi(t ■ Ti))<T ■ n cm (Tj • fl c 

+ (C + Cl(T • 7i))£T + (C + C*i(r ■ Ti))ai 

+ (E + Ex(r ■ T i ))a- ■ qcmo-i ■ q cm 

+ (F + Fi(r ■ Tj))(j ■ QcmO-j • Q cm , (4.5) 

W — W k y W 

Qcm 7j~} j~ Ti "cm 7j~ T~} f ) vicm — Qcm X n cm , l^.OJ 

l^cm — K cm| |K C m x k cm\ 

and k cm (k£ m ) is the projectile (ejectile) momentum in the NN cm. system. 

The relation between the NN f-niatrices and the NN scattering amplitudes in the NN cm. frame is 



where 



^NN^lab, 9cm) = ^/nN (-f^lab, 9cm), (4.7) 
-47T 

v/M 2 + fc, 2 



— 47T 

V= * ,„ • (4.8) 



cm 



Therefore the in (3.9) are given by 



t U u — v(Ao + Ai(tq ■ Ti)), (4.9) 

tnn=V(Bo+B 1 (T Q -T l )), (4.10) 

tm. =^(Co+C 1 (r -r l )), (4.11) 

tun=v(Cv+C 1 {T -T i j), (4.12) 

t 9g = 7?(£;o+£i(ro-T i )), (4.13) 

V = 7 7( i? o+^i(To-r,)), (4.14) 

and all others are 0. 

Calculating the expectation values of the isospin operators, we get the amplitudes for proton-proton 
scattering as 

t ua = r)(A +A 1 ),---. (4.15) 

For neutron-neutron scattering these are the same as for proton-proton scattering. 

Other types of NN collisions are proton-neutron charge exchange collisions and proton-neutron non- 
charge exchange collisions. Their amplitudes are 

t uu = 2r)Ax,---, (4.16) 

t uu = r)(A -A 1 ),---, (4.17) 

respectively. 

We note that the isospin operators of the struck nucleons are included in the response functions. 
Hence we use the i-matrix amplitudes for the charge exchange collision as 

t uu = V2rjAi,---, (4.18) 



instead of ( 4.16 ) 
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4.4 Response functions 



We obtained a simple expression for the response function in (2.31). We take account of the weight 
function cb(u>) through a complex potential for the single particles, which is of a Woods-Saxon-type 
central potential with the derivative form spin-orbit potential as 



U'(r) = -(V + iW)- 



1 



1 + exp 



r-Ro 



V la 



a 

fr~R so \ 
exp 

V a so / 



run J aso 



1 + exp 



r - Rso 



r(l • s) + Vboui, 



(4.19) 



where the radius parameters are r = rgo = rc ou i — 1-27 fm, the diffusenesses arc a = ago = 0.67 fm, 
the potential depths are determined by the binding energy of the last occupied single particle level, and 
the imaginary part of the potential is W — 5.0 MeV. J25| The depth of the 1 • s part is chosen to be 6.5 
MeV (10.0 MeV) for 12 C ( 40 Ca). @ The response functions used in the DWIA calculation depicted 
in Figs. 1 and 2 include the effects of the RPA correlation with the A degree of freedom as well as the 
effective mass and the spreading width of hole states. However, here we calculate the response functions 
with no correlations in order to extract the effect of the nuclear reaction mechanism exclusively. 

We extended the response functions to nondiagonal form with respect to the momentum transfer q 
in the two-step formalism. The spin scalar isoscalar response function in the nondiagonal form is defined 
as M 



B=£0 













i 




3 



$ ) 5{u mt - E 1 ^) 



(4.20) 



The response functions which include the spin operators are generally written in the form 
iWq,qV nt ) 





ue 



$o S(co im - Eg*). 



(4.21) 



Here the response function has also been extended to the nondiagonal form with respect to the spin 
direction. 

If the spin-orbit force is negligible and the target nucleus is spin-saturated, the single particle states 
can be written as \n, l,m,s = 1/2, //) with the principal quantum number n, the angular momentum I, its 
third component m, the spin s, and its third component /i. Then the response function can be separated 
into the spin part and the orbital part as 

i?u'u(q,qV nt ) ee X) K/*fck-u'|A«p)(A*pk-u|/*fc)] 

n p l p m p n h l h m h (IpHh 

x [(n ft / ft m ft |e lq '' r |np^mp)(n p /pm p |e" 4q ' r |n/ l // l m /l ) 6(u iut - E^f), 

(4.22) 

where p and h denote a particle and a hole, respectively. The spin part becomes Tr[<r • u'er • u] = 2u' • u, 
while that of the response function without the spin operators becomes the trace of the unit matrix. 
Thus we get the relation 

i?a<a(q, q'; cu iat ) = (u' • u)i?(q, q'; W int ). (4.23) 



This relation is plotted in Fig. 5, and we judge it reasonable. 
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If only one spin operator is included in the response function, then the spin part becomes Tr[<7-u] = 0, 
and the response functions become 



i? fi (q,qV nt ) = 





$o } 6(u int - Eg) 







Ub 




i 






3 



$o > 6(u> 



int /pint > 
Eg ) 



(4.24) 

We cannot calculate the response functions for charge exchange collisions with these relations. There- 
fore we must prepare two more response functions. One is for the (p, n) collision and the other is for the 
(n,p) collision. The difference among these three response functions is seen in Fig. 6. These response 
functions also satisfy the relations 

(4.25) 
(4.26) 



where 



i?+(q, q';w int ) 



iT (q,qV nt ) 




(q,qV nt ) = 


= (u' 


■u)R+ 


(q,q';^ int ), 


(q,qV nt ) = 


= (u' 


■ u)R- 


(q,q';^ int ), 


i 




) ($B 


3 






) ($B 




i 









(4.27) 

with r± = (t x ± iTy)l\pl. The response functions i? + and i?~ are for the (p,n) and (n,p) collisions, 
respectively. With these three response functions, it is sufficient to calculate the various response functions 
that appear in ( [4.3[ ). Using the 'spectator assumption', we assume that the target is in the ground state 
of 12 C ( 40 Ca), even when (n,n') or (n,p) collision occurs in the second step. 

We must consider the struck nucleons when w e use the response function for (p,p') collisions. If we 
write the response function for (p,p') collisions as ( 4.20|) , the struck nucleon can be a proton or a neutron. 
However, when we calculate the cross sections, we distinguish a proton from a neutron, as shown in Figs. 3 
and 4. Therefore we must divide the spin scalar isoscalar response function by 2 when we apply this 
response function to (4.3). 



4.5 Cross sections 



Using the spin-dependent form of the cross sections in (3.14), the one-step unpolarized, the spin longitu- 
dinal and transverse cross sections are written as 



/ = i 4rElW(q)l 2 ^(qX nt )> 



ID q = K ^-\t qi} (c l )\ 2 R qil ( q. .•""). 
/£ p = i^Mq)| 2 i? pp ( q y nt ) 

MR 



(4.28) 

(4.29) 
(4.30) 



respectively. Here we have used the approximation q ~ q cm which was justified in Rcf. |lj. The 
summations of fi and /2 run over only (ji, p,) = (u,u), (n,n), (n, u) , (u, n ), (q,q) and (p,p), and we have 
used the relation i? p ^(q, u lnt ) = R(q, uj lnt ), which is derived from (4.23). 
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We now explain the method of calculating the two-step cross sections given in ( 3.19| )-( |3.21 ). The 



(2) 

quantity X;,' , has four NN i-matrices, each of which generates 6 terms. Therefore we encounter 



6 (= 1296) terms for each isospin-dependent process explained in §4.1. However, with the approximation 



(4.24) we can eliminate many terms, and finally we have to take only 400 terms (see Fig. 7). We see 



that the response functions are common except the scalar products that appear from the spin parts (see 



( [4.23D and so on). 

Next we carry out the summation of 400 products of four NN i-matrices with the scalar products. 
The traces of the spin operators of the incident nucleon are easily calculated as, for instance, 

-Tr(<7 9 cro M2 cTo Ml ) = s • q, (4.31) 

iTr(crp<7 0M2 cr 0Ml ) = s • p, (4.32) 

with the relation 

o"o M2 (To Ml = I sq + cr • s, (4.33) 

where Sq is a real number and s is a complex vector. After the summation we calculate the five-fold 
integration of the sum multiplied by the spin scalar response functions. 

In fact, the spin scalar response functions are also common if charge exchange collision occurs in the 
same step in the (p, n) reactions. As for the (p,p') scattering, they are common to the four processes that 
have no charge exchange collision (see Figs. 3 and 4). Hence, in a practical calculation, we also carry 
out the summation of the products of four NN i-matrices with the scalar products in these processes, in 
which the spin scalar response functions are common, and after the summation we calculate the five-fold 
integration. 



5 Results 

In Fig. 8 we display the double differential cross sections of the (p,p') scattering and the spin longitudinal 
and the spin transverse 12 C(p,n) reactions at 494 MeV as functions of the energy transfer in the 
laboratory frame. In this calculation the momentum transfer q mt in the intrinsic frame of the target is 
1.55 (1/fm) (qiab in the laboratory frame is 1.69-1.72 (1/fm)). The scattering angle corresponding to this 
momentum transfer is about 18°. (ll], |lj, 

In the 12 C(p,p') scattering, the peak of the one-step cross section is around uj\ a b — 60 MeV in this 
no-correlation calculation, and the contribution of the two-step processes there is only about 0.5% of the 
one-step process. The two-step cross section is still increasing around o>i a b = 125 MeV, and it becomes 
about 6.0% of the one-step process there. These results can be seen in panel (d), which indicates the 
ratios of the two-step cross sections to the one-step cross sections. The contribution of the two-step 
processes to the 12 C(p,p') scattering is negligibly small. However, in the 12 C(p, n) reactions, there are 
considerable contributions from the two-step processes. In the spin longitudinal 12 C(p, n) reaction, the 
contribution of the two-step processes is 3.5% of the one-step process at oj\ a b = 60 MeV and 28% at 
cji a b = 125 MeV. For the spin transverse 12 C(p, n) reaction, the two-step contribution becomes larger 
than that for the spin longitudinal reaction, and at oj\ a b = 60 and 125 MeV the ratios become 9.0% and 
68%, respectively. From these results we can say that the two-step processes become more effective as 
the energy transfer increases in the quasi-elastic region. 

Figure 9 displays the double differential cross section at 346 MeV of the same reactions. The mo- 
mentum transfer q lnt is 1.52 (1/fm) (qi a h is 1.66-1.68 (1/fm)). It corresponds to a scattering angle of 
22°. |14j The contributions from the two-step processes become larger than thoseat 494 MeV. In the 
(p,p') scattering, the ratio at wi a b = 60 MeV is less than 1% and at u>\ a b — 125 MeV it is about 11%. In 
the spin longitudinal 12 C(p, n) reaction, the ratio is close to 40% at ^i a b = 125 MeV, while that of the 
spin transverse 12 C(p, n) reaction exceeds 75% there. The contribution to the spin transverse reaction is 
much larger than that to the spin longitudinal reaction, as in the case of 494 MeV. 

We assume that the i5-term and the F-term in the NN i-matrices mainly determine the spin lon- 
gitudinal and the spin transverse cross sections, respectively, in the two-step processes. The difference 
between the two-step contributions for the spin longitudinal and the spin transverse cross sections may 
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be due to the momentum transfer dependence of the E-term and the F-term which are shown in Fig. 10. 
The dominant real part of the E-term changes sign at 0.7 (1/fm), while the F-term never becomes 0. The 
square of NN i-matrices for proton-proton and non-charge exchange neutron-proton scattering decreases 
monotonically as a function of the momentum transfer. Thus the contribution from the F-term becomes 
smaller than that from the F-term through the integration of the momentum transfers. 

The results for 40 Ca at 494 and 346 MeV are given in Figs. 11 and 12, respectively. Since 40 Ca 
is a spin-saturated nucleus, the relation (4.23) becomes a better approximation. We guess that the 
results are more reliable than those for 12 C. The two-step contributions in comparison with the one-step 
contributions are larger than those for the 12 C case. The momentum transfer q lnt for the reactions at 494 
MeV is 1.655 (1/fm) (q\ a h is 1.69-1.71 (1/fm)). For the (p,p') scattering and the spin longitudinal and 
spin transverse 40 Ca(p,n) reactions, the ratios are 7.3%, 45% and 87% at cj lab = 120 MeV, respectively. 

At 346 MeV the momentum transfer q mt is 1.625 (1/fm) (qi a b is 1.66-1.68 (1/fm)). The ratios reach 
13% for the (p,p') scattering, 49% for the spin longitudinal (p, n) reaction and 88% for the spin transverse 
reaction at cj[ a b = 120 MeV. We can say that the contribution from the two-step processes is larger at 346 
MeV than at 494 MeV only for the (p,p') scattering in the 40 Ca case. Here we have confirmed again that 
the two-step contribution to the spin transverse reaction is much larger than that to the spin longitudinal 
reaction. 



6 Summary and conclusion 

We constructed a formalism for the two-step direct reaction within the framework of the plane wave 
approximation in order to compare its contribution with the one-step process. With this approximation 
we factored the nucleon-nucleus T-matrix into the NN t-matrix and the transition density in the one-step 
formalism. In the two-step formalism we expressed the two-step cross section with a nondiagonal response 
function with respect to the momentum transfer q. 

In the Green's function we removed the effect of absorption, because the motion of the projectile and 
the ejectile was described with plane waves. We applied the on-energy shell approximation to the Green's 
function. 

With this formalism, we calculated the cross sections of the 12 C, 40 C&(p,p') scattering and the spin 
longitudinal and spin transverse cross sections of the 12 C, 40 Ca(j3, n) reactions at 346 and 494 MeV. The 
scattering angles were set to 22° and 18°, respectively, for comparison with the experiments at LAMPFp3| 
and RCNP.f4| 

In (p,p') scattering, the contributions from the two-step processes were found to be small. However, 
there are appreciable contributions from two-step processes in the (p, n) reactions. In the spin longitudinal 
(p, n) reaction, the two-step contribution is about 30%-50% of the one-step cross section, while the ratio 
for the spin transverse (p, n) cross sections become about 70%-90% around wi a b = 120 MeV. We believe 
that the difference between the two-step contributions for the spin longitudinal cross section and the 
spin transverse cross section is due to the difference between the momentun transfer dependences of the 
F-term and the F-term in the NN scattering amplitude. 

We found that the contributions of the two-step processes become larger as the energy transfer in- 
creases in the quasi-elastic region. We also found that the two-step contributions in comparison with the 
one-step contributions are larger in the 40 Ca case than in the 12 C case. 

It has been reported that in DWIA calculations the cross sections of the spin longitudinal (p, n) 
reactions are underestimated beyond the quasi-elastic peak region, and for the spin transverse (p, n) 
reactions, the calculations amount to only half of the experimental results. We give the DWIA results 
multiplied by the ratios of the sum of the one- and the two-step cross sections to the one-step cross section 
for 12 C in Fig. 13 and for 40 Ca in Fig. 14. From these figures we can clearly see that the theoretical 
results including the two-step contribution are closer to the experimental results than the DWIA results, 
particularly for ID q at 494 MeV. 

However, our results are still insufficient for explaining all of the discrepancy. In those regions where 
the cross sections are still underestimated, the response functions that include 2p-2h correlations may have 
some effect. Further multi-step processes must be investigated in these regions. A two-step calculation 
with full distortion is of course needed to obtain a more quantitatively reliable conclusion. 
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Appendix 

We used the method of Love and Franey j22| to calculate the NN i-matrix. It is written in the form 

tMKub,q) = [Vsofa) - V S C (Q)]P S0 + [V s c E (q) + V&(Q)]Pse 
+ [Vroto - Vto(Q)]^to + [Vf E (q) + K r c E (Q)]P TE 

+ l -[QV^°{q) + qV LSO (Q)](°0 +<7i) • hP TO 

+ l -[QV^(q) - qV LSE (Q)](a + a, t ) ■ nP TE 

- [V™°(q)S?™(q) - V™°(Q)S?™(Q)} 

- [V™(q)Sf™(q) + V™(Q)Sf?™(Q)}, (A.l) 



where 



and 



q = 2fe cm sin(0/2), (A.2) 
Q = 2fc cm cos(#/2), (A.3) 



(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 



Pso 


= Ps= 


=o-P*=o — 


1 - c 


■CTil — 


t ■ n 


4 




4 


PsE 


= P S = 


-.vPt=i = 


1 - c 


■ <?i 3 + 


T ■ Ti 


4 




4 


Pto 


= P S = 


=liw = 


3 + cr 


• <Ti 3 + 


t ■ n 


4 




4 


Pte 


= Ps= 


:1-F*=0 = 


3 + (T 


• Uj 1 - 


TO • Tj 


4 




4 


S?? D (q) 


— S12 


(q)Pt=i -- 


= (3er ' 


go* • g 


- (T ■ 


S??™(q) 


— S\2 


(q)Pt=o -- 


= (3cr • 


• q 


- To ■ 



• T • Tj 



4 

i - t • n 



4 

with the projectile wave number k cul in the NN cm. system. The potential depths are given as 

^-^ CT 1 <A ' U) 

where fe = q or Q, s is the spin singlet (S) or spin triplet (T), and tt is the parity in the NN system. The 
summations over i are taken over several ranges. It is these values, i3 V A i LS7r and I^ TN7r that are given 
in the table in Ref. p3[ . 
Defining 



r so 




-^so( 




= Vg(q) - 


+ vg( 


E TO 




v TO 
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^SO= l [Qt? LSO (g) + gt? LSO (Q)]j 

t LSB s l [ g^LSE ((?) _ g ^LSE (Q)]) 
t TNO = ^0 (g)j ^TNO (Q)j 

t TNE Ey™( g ),y TNE (Q), (A.13) 



and noting Eqs. ( |4.5| ) and (4.7), one obtains the relations 

rjA = i(ig Q + 3t£ E + 9t? + 3t? E ), (A.14) 

vAi = -^(-*so + *se + 3i?o " 3t? E ), (A.15) 
rjBo = *so — 3<g E + 3txo + *te); 

+ I(3t TNO (g)+t TNE (g)-3t TNO (g) + i TNE (Q)), (A.16) 

^1 = Jq^SO ~ ^SE + *TO — *Te) 

+ \(t TN °(q) t™(q) - i TNO (Q) - t TNE (Q)), (A.17) 

? ? C = i(3t LSO +t LSE ), (A. 18) 

??Cl = i(t LSO -i LSE ), (A.19) 

VE = VB + \(-9t TNO (q) - 3i TNE (g)), (A.20) 

VEi = vBi + J(-3t TNO (g) + 3i TNE (g)), (A.21) 

vFo = vBo + J(9t TNO (Q) - 3i TNE (Q)), (A.22) 

TfFx = rjBt + i(3t TNO (Q) + 3i TNE (Q)). (A.23) 

The corresponding t-matrix in the nucleon-nucleus system can, to a good approximation, be written 
by replacing Q with the wave number of the projectile in the nucleon-nucleus system ki, which is expressed 
as 

A + a 



k( = M Z A(3 



1 + 13 



«=^r> P^nrr^^ (A. 24) 



if lab B= 4aA 
2M ' P ~ {A + 1) 

where A is the mass number of the target nucleus. 

Moreover, for calculating nucleon-nucleus scattering, an A dependent kinematic modification is re- 
quired. This modification is provided by the transformation of the t-matrices given by 

tNA = - 2 — £NN) (A. 25) 



' -M 2 (l + a), ef = M 2 + kf, e 2 = M ' 2 + ( T ) > ( A - 26 ) 



cm 



where £i (et) is the total energy of the incident (target) nucleon in the nucleon-nucleus system, and e c 
is the total energy of the incident nucleon in the NN system. 
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Figure 1: The spin longitudinal cross section ID q and the spin transverse cross section ID p of the 
12 C(p, n) reaction at the incident energy 494 MeV calculated in DWIA. [17] Experimental data are taken 
from Ref. [13]. 
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Figure 2: Same as Fig. 1 but at 346 MeV.[14] 
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Figure 3: Four processes of the two-step (p, n) reactions. The motion of the projectile and the ejectile 
is drawn on the left side of each process. The particles drawn on the right side depict the motion of the 
struck nucleons in the first-step and the second-step collisions. 
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Figure 5: Examination of (4.23) in the case that u = q and u' = q'. The solid line indicates the 
right-hand side of (4.23), and the dashed line indicates its left-hand side. The momentum transfers 
are q = q' = 1.6 (1/fm), = 0' = 30° for the left panel and = 0'= 90° for the right panel, and 
<P = 0°, 4>' = 180°. 




Figure 6: Spin scalar response functions of the diagonal part with respect to the momentum transfer 
q (a) for 12 C and (b) for 40 Ca. The solid line indicates the response function for (p,p') scattering, the 
dashed line indicates that for (n,p) reactions, and the dotted line indicates that for (p, n) reactions. The 
correlation is not included. 
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Figure 7: The combinations of four amplitudes in NN i-matrices that do not disappear. The two ovals on 
the left side of each panel represent the NN i-matrices in the first and the second step. Those on the right 
side represent their complex conjugates. The lower box is the response function in the first step, and the 
upper one is that in the second step. The letters A, B, C, E and F are the amplitudes in (4.5). If t(q\) 
is the A-tevva or the C-term, which do not include the spin operators for struck nucleons, then tf(q[) also 
has to be the A-term or the C-term. Otherwise, the response function in the first step includes only one 
spin operator and it becomes 0. From the four panels 16, 64, 64 and 256 terms appear, respectively, and 
in all 400 terms appear. 
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(b) 12 C(p,n) ID t 
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Figure 8: (a) Double differential cross sections of the 12 C(p,p') scattering, (b) The 12 C(p, n) spin lon- 
gitudinal reaction at 494 MeV. (c) The 12 C(p, n) spin transverse reaction at 494 MeV. The dashed and 
dotted lines indicate the one-step and the two-step cross sections, respectively. The solid line is the sum 
of these two cross sections. The momentum transfer gj a b in the laboratory frame is about 1.70 (1/fm). 
(d) The ratios of the two-step cross sections to the one-step cross sections. The solid line, the dashed 
line, and the dotted line indicate the (p,p') scattering, the spin longitudinal (p, n) reaction, and the spin 
transverse (p, n) reaction, respectively. 
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Figure 9: Same as Fig. 8, but at 346 MeV. The momentum transfer gi a b is about 1.67 (1/fm). 
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Figure 10: The isovector spin longitudinal and spin transverse NN t- matrices. 
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Figure 13: The DWIA calculations multiplied by the ratios for the 12 C(p, n) reactions. The left panels are 
the cross sections for the spin longitudinal reactions, and the right panels are those for the spin transverse 
reactions. The dashed lines indicate the DWIA results, and the solid lines indicate those multiplied by 
the ratios. 
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Figure 14: Same as Fig. 13, but for 40 Ca 
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